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Abstract. We study quantum transport in honeycomb lattice ribbons with either armchair or zigzag edges. 
The ribbons are coupled to semi-infinite linear chains serving as the input and output leads and we use 
a tight-binding Hamiltonian with nearest-neighbor hops. For narrow ribbons we find transmission gaps 
for both types of edges. The center of the gap is at the middle of the band in ribbons with armchair 
edges. This symmetry is due to a property satisfied by the matrices in the resulting linear problem. In 
ribbons with zigzag edges the gap center is displaced to the right of the middle of the band. We also 
find transmission oscillations and resonances within the transmitting region of the band for both types 
of edges. Extending the length of a ribbon does not affect the width of the transmission gap, as long as 
the ribbon's length is longer than a critical value when the gap can form. Increasing the width of the 
ribbon, however, changes the width of the gap. In armchair edges the gap is not well-defined because of the 
appearance of transmission resonances while in zigzag edges the gap width systematically shrinks as the 
width of the ribbon is increased. We also find only evanescent waves within the gap and both evanescent 
and propagating waves in the transmitting regions. 



PACS. 73.23.-b Electronic transport in mesoscopic systems 
materials and structures - 05.60.Gg Quantum transport 



73.63.-b Electronic transport in nanoscale 



1 Introduction 

Graphene is a single layer of carbon atoms in a honeycomb 
lattice and serves as the basic building block of fullerenes 
like carbon nanotubes pQ. It has recently been found to 
be stable under ambient conditions [5J[3] an d has unique 
charge carrier properties ifi^&i. Graphene nanoribbons 
are quasi-one-dimensional honeycomb lattice sheets with 
finite widths. Because of the width termination charge car- 
riers become confined and the appearance of an energy 
gap is expected f7tf8|[9|fl^ | flT 1 [l ^ [l3 l ll4 p 5j . For a recent re- 
view see reference |16j . Graphene nanoribbons have been 
fabricated using lithography [17], lithography and etch- 
ing [18], epitaxial growth and lithography [19], chemical 
growth [20], and deposition and etching [21] techniques. 
Measurements on these nanoribbons found a conductance 
gap that scales inversely with the width of the ribbons (TTJ, 

in. 

Possible edge types of regular graphene nanoribbons 
can be armchair or zigzag edges, on both sides of the rib- 
bon. Theoretical studies using tight-binding band calcula- 
tions based on the 7r-states of carbon [5|f5I[TTJ|[TT] and stud- 
ies on the solution of the two-dimensional free massless 
Dirac equation using specific boundary conditions 12,13, 
1 14] found different transport properties for ribbons with 
armchair edges than those with zigzag edges. Ribbons 



with armchair edges can be either metallic or semicon- 
ducting depending on the width of the ribbon while rib- 
bons with zigzag edges are metallic. Furthermore, the ex- 
istence of special edge states localized at the zigzag edges 
gives rise to flat bands at the Fermi level in ribbons with 
zigzag edges. On the other hand, ab initio calculations 
on graphene nanoribbons with hydrogen passivated edges 
found band gaps in both armchair and zigzag edges [7]. 
Another density functional theory calculation, however, 
found that hydrogen termination can affect the properties 
of the ribbon [15] . 

In molecular electronics the search for conducting molec- 
ular wires that can be used in the construction of devices 
such as field effect transistors has been actively pursued. 
Wires made of linked thiophene units 22,23,24,25 and 
alternating thiophene and ethynyl units [26] have been 
synthesized with lengths up to 11 nm. The synthesis of a 
chain of copper ions surrounded by helical ligand strands 
has also recently been reported [27] . Such molecular wires 
can be coupled to electrodes such as cut nanotubes [28, 
[29] or gold break junctions [30ll31j. 

In this paper we study quantum transport in hon- 
eycomb lattice ribbons using a tight-binding model with 
nearest-neighbor hops. The ribbons have either armchair 
or zigzag edges. In contrast to previous theoretical and 
numerical studies, however, we couple the ribbons to semi- 
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infinite linear chains that serve as the input and output 
leads. Such a system may be realized experimentally as a 
nanoribbon of graphene in contact with molecular wires 
serving as the leads. We vary both the length and the 
width of the ribbons. For a chain of hexagons with either 
armchair or zigzag edges and coupled to input and output 
leads we find the existence of a transmission gap, as long as 
the ribbon is longer than a critical length below which the 
gap does not form. The center of this gap coincides with 
the middle of the band in ribbons with armchair edges 
but not with zigzag edges. This symmetry arises from a 
property of the matrices in the resulting linear problem. 
Extending the length of the ribbon does not affect the 
width of the transmission gap. Extending the width of 
the ribbon, however, affects the width of the gap. In rib- 
bons with armchair edges the appearance of transmission 
resonances obscures the exact width of the gap while in 
ribbons with zigzag edges the gap systematically shrinks 
as the width of the ribbon is extended. 



2 The model 

We model the transport of a quantum particle using the 
tight-binding Hamiltonian 

N„ 

n = J2 e * \<Pi) (<Pi\ + (fj\ + Ifj) (<Pil}> (!) 

i=l <ij> 

where the \<Pi)'s are tight- binding basis functions centered 
on site i, N s is the total number of lattice sites available for 
the particle to hop into, and the sum in the second term 
is only over nearest-neighbor sites. Randomly varying the 
on-site energies e$ identifies the model as the Anderson 
model of localization [32] • Alternatively, setting the ej to 
a constant, choosing Vij to be a non-zero constant only for 
nearest neighbor sites, and embedding the transport in 
a disordered site-percolation cluster identifies the model 
as the quantum percolation model 33,34]. Both the An- 
derson model and quantum percolation involve disorder. 
In this paper we follow the quantum percolation model, 
i.e., we set = and let Vij be a non-zero constant only 
for nearest neighbor sites. The model, however, does not 
involve disorder. Particle transport is through perfectly 
ordered honeycomb lattice ribbons. 
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Fig. 1. (color online) Hexagons in (a) armchair and (b) zigzag 
orientations. Semi-infinite input and output linear chains are 
connected to the hexagons. Beside each site is its label. Sites 
labelled and JV r + l are the input and output supersites, re- 
spectively. Ni is the number of sites within the hexagon. 



Shown in Fig. [T] are hexagons oriented so that they 
form either armchair or zigzag edges in a honeycomb lat- 
tice ribbon. In a ribbon containing N r sites, the site with 
label in the input chain and label N t + 1 in the output 
chain are the only sites that are directly connected to the 
ribbon. We call these sites the input and output super- 
sites, respectively. The connections between these super- 
sites and the ribbon can be single-channel as shown in 
Fig. Ufa) or multi-channel as shown in Fig.QJb). The dis- 
tance between nearest-neighbor sites is set to ao = l, thus 
setting the length scale of the model (although the length 
scale is actually inconsequential since Ti has no explicit 
length dependence). 



N = 1, m = 3, armchair N = 1, m = 3, zigzag 




Fig. 2. (color online) Honeycomb lattice ribbons of various 
widths. N is the width label of a ribbon and m is the number 
of hexagons in one hexagonal chain. Connections between the 
ribbon and the linear chain leads are also shown. 

Shown in Fig. [5] are honeycomb lattice ribbons with 
either armchair or zigzag edges and with varying widths. 
The N = 1 ribbon is a chain of hexagons coupled to the 
input and output leads, m counts the number of hexagons 
in one hexagonal chain. The N=2 ribbon consists of two 
N = 1 ribbons placed side by side. N thus directly de- 
termines the ribbon width while m determines the ribbon 
length. 

We have used semi-infinite linear chains as the input 
and output leads to mimic linear molecular wires. Because 
of this linearity, an incident particle having a specific en- 
ergy can have only one mode of either a propagating or 
evanescent wave 35 impinging on the ribbon. We thus 
determine the effects of small honeycomb lattice ribbons 
of varying widths and lengths on the transmission of such 
single-mode waves. Recent studies have been done mod- 
elling the leads as normal metals having an aggregate 
number of modes incident on large honeycomb lattice rib- 
bons [36 , 37, 38J . The leads were modelled as either square 
or honeycomb lattices. It was found that for sufficiently 
large ribbons the geometrical details of the leads were not 
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critical [36]. At the Dirac point of the ribbon evanescent 
waves are dominant while away from this point propagat- 
ing waves dominate the transport |37j . In this study, in 
contrast, we deal with leads having only one mode inci- 
dent on a small honeycomb lattice ribbon and see that the 
transmission of this incident wave depends strongly on the 
width, length, and edge type of the ribbon. 

To determine the hopping transport properties of par- 
ticles through honeycomb lattice ribbons the Hamiltonian 
in Eq. (JlJ is cast in a matrix representation with the \<fi)'s 
as the basis. We write the eigenvector as which is 
a column vector with elements ipi as the coefficients of 
the chosen basis. Because the linear chain leads are infi- 
nite the equivalent matrix problem is infinite-dimensional. 
The contribution of the leads (except for the supersites), 
however, can be factored out since each site in the chain 
only interacts with its nearest-neighbor sites to the left 
and right, thereby resulting in predictable tridiagonal el- 
ements in the Hamiltonian matrix. We follow the method 
proposed by Daboul et al. [39j by using the scattering 
boundary condition ansatz 



e mq +re- mg , 

t e i(n-N r -l)q 



left input lead, 
right output lead, 



(2) 



where r and t are the reflection and transmission ampli- 
tudes, respectively. Note that along the left input lead 
n<0 while along the right output lead n> N r . The ansatz 
states that an incoming plane wave and a reflected wave 
with amplitude r are in the input lead and a transmitted 
wave with amplitude t is in the output lead. An alterna- 
tive but equivalent approach is to use the non-equilibrium 
Green's function technique to determine the transmission 
00]. 

In general the hopping amplitudes along the left and 
right leads and the ribbon can be different. In this paper 
we consider the left and right leads to be the same and 
thus, the leads have the same hopping amplitudes t>i e ft = 
bright = ^icad- Let v be the hopping amplitude within the 
ribbon. 

Using the ansatz we can reduce the original infinite- 
dimensional problem into a finite one. Along both leads 
we get the dispersion relation 



E 



(3) 



Choosing a numerical value for E and v\ cat \ sets the value 
of q. Also, because ofEq.[3l the value oiE/vi ea d is bounded 
within the interval [—2,2]. For the ribbon, the problem 
reduces to a finite linear problem of the form Aip = b. In 
the armchair hexagon shown in Fig.QJa), for example, we 
get 

/-a+i/3 wioad \ 
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(4) 



where a = E/2, (3 = y/^ cad -E 2 /2. 

We can determine tp exactly by carefully taking the 
inverse of A and then multiplying it to b. To do this we 
numerically decompose A using singular-value decompo- 
sition so that A = UlVt, where U and V are unitary 
matrices and Z is a diagonal matrix [31]. The inverse can 
then be determined using the properties of unitary and 
diagonal matrices. Once we determine if) we can then cal- 
culate the transmission coefficient by T — \ t\ 2 . This numer- 
ically exact approach has also been employed in the study 
of localization in disordered clusters in two-dimensional 
quantum percolation [42U43j . 

The connections between the input and output super- 
sites and the ribbon are reflected in the outermost rows 
and columns of A. The inner rows and columns reflect the 
connections between sites within the ribbon. And so for 
the zigzag hexagon shown in Fig. [JJb) the corresponding 
matrix A will have almost the same central entries as that 
shown in Eq. ((4]) except for the extra entries along the 
outer rows and columns because of the additional chan- 
nels available between the supersites and the hexagon. 

A merit of determining all elements of ip exactly is that 
the components tpt of the wave function at all of the sites 
in the ribbon are also determined. In Sec. [3] we show plots 
of the wave function density in N — l, m = 8 ribbons with 
armchair or zigzag edges. 



3 Numerical results 

Depending on the atomic composition of the ribbon and 
the leads, v and wi ca d are in general not the same. In our 
simulations, however, we set the values of v and wi ea d to 
one. We vary the energy from E = —2 to E = 2 in energy 
steps of AE = 0.01. For each value of E the complex- 
valued matrix A is constructed by noting the connections 
between the sites, including the supersites. The function 
ip is then numerically determined via the singular-value 
decomposition of A. The transmission coefficient T is then 
calculated from \t\ 2 . 

Shown in Fig.[3]are plots of the transmission coefficient 
T as a function of the energy E of the particle m N — l rib- 
bons. Ribbons with armchair edges are shown in Fig. [3Ja) 
while those with zigzag edges are shown in Fig. [3fb). The 
length of the ribbons are m = 2 and m = 8 hexagons. 
For the ribbons with armchair edges we see symmetry be- 
tween the positive and negative energy regions, i.e., with 
respect to the middle of the band. This symmetry in the 
sign of the energy arises from a condition satisfied by the 
A matrix for the ribbon with armchair edges and is dis- 
cussed in Sec. |3J We also see two features as the length of 
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the ribbons are increased from m~2 (red dots) to m = 8 
(black dots): (a) the number of transmission oscillations 
increases proportionally, and (b) a transmission gap devel- 
ops. The transmission oscillations are reminiscent of the 
transmission resonances found in the square lattice case 
[44] . The transmission gap develops as the length of the 
ribbon is extended. Although T has a non-zero minimum 
at the middle of the band for m = 2, this minimum falls 
to zero transmission for the first time when m = 6. Ex- 
tending the ribbon length further to m — 100 increases 
the number of transmission oscillations and sharpens, but 
does not change the width, of the gap. 




Fig. 3. (color online) Transmission coefficient T as a function 
of the particle's energy in honeycomb lattice ribbons of width 
TV — 1. (a) Armchair edges with m = 2 for filled (red online) 
circles and m = 8 for open circles. The lines drawn are not 
analytical fits but simply connect adjacent datapoints to aid 
the eye. (b) Zigzag edges with m = 2 for filled (blue online) 
circles and m — 8 for open circles. 

For the ribbons with zigzag edges we see in Fig. [3](b) 
that the number of transmission oscillations grows and a 
transmission gap develops as the length of the ribbon is 
extended from m = 2 (blue dots) to m = 8 (black dots), 
similar to the case for armchair edges. However, the cen- 
ter of the transmission gap is not at the middle of the 
band, which is now slightly transmitting. And in contrast 
to the armchair edges case, the number of oscillations in 
the E > region is half those in the E < region. The ab- 
sence of symmetry with respect to the middle of the band 
in the transmission for ribbons with zigzag edges can be 
attributed to the extra connections between the input and 
output supersites and the ribbon, as discussed in Sec. 3J 




Energy Energy 

Fig. 4. Transmission coefficient T as a function of the particle's 
energy E in honeycomb lattice ribbons with armchair edges. 
The ribbons have length m = 100. The widths are (a) TV = 2, 
(b) TV = 3, (c) TV = 4, and (d) TV = 5. The transmission gap is 
symmetric with respect to the middle of the band. 




-2-1012-2-1012 
Energy Energy 



Fig. 5. Transmission coefficient T as a function of the particle's 
energy E in honeycomb lattice ribbons with zigzag edges. The 
ribbons have length m = 100. The widths are (a) TV = 2, (b) 
TV = 3, (c) TV — 4, and (d) TV = 5. The transmission gap shrinks 
systematically as the ribbon's width is increased. 

Shown in Figs. 2] and O are plots on how the transmis- 
sion coefficient behaves as a function of the energy of the 
particle as the width of the ribbons is varied from TV = 2 
to TV = 5. The ribbons have length m= 100 and have arm- 
chair edges in Fig. [4] and zigzag edges in Fig. [5] For the 
ribbons with armchair edges the transmission gap is sym- 
metric with respect to the middle of the band. In contrast, 
for the ribbons with zigzag edges the center of the trans- 
misison gap does not coincide with the middle of the band. 
In addition, the width of the gap decreases systematically 
as the width of the ribbon with zigzag edges is increased. 
This is not the case for the ribbons with armchair edges. 
Taking a closer look at Fig. BJb) we see that transmission 
resonances appear in the region where there is a gap in 
the narrower ribbon shown in Fig. HJa). As the ribbon's 
width is increased transmission resonances appear within 
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Table 1. Values of the fitting parameters found for the data 
shown in Fig. [6] for ribbons with zigzag edges. 



linear fit: 

A g = (-0.082 ± 0.010) N + (0.917 ± 0.065), \R\ 2 = 0.928 
power- law fit: 

A g = (1.061 ± 0.058) at(-o-60o±o.o49) ^ = q ggy 
exponential fit: 

A a = (1.175 ± 0.038) e (- 0197±0 009 » N , \R\ 2 = 0.992 



the gap and thus obscures the determination of the exact 
width of the gap. 

We thus see in narrow ribbons with either armchair or 
zigzag edges a transmission gap dividing the transmission 
curve into two regions. As the width of the ribbon is in- 
creased the gap in ribbons with armchair edges is obscured 
by the appearance of transmission resonances within it 
while the gap in ribbons with zigzag edges systematically 
shrinks. 




1 1 1 1 1 1 

2 4 6 8 10 

N 

Fig. 6. (color online) Transmission gap, A g , as a function of 
the width N for m = 100 ribbons with zigzag edges. Linear 
(red), power-law (green) and exponential (blue) fitting curves 
are also shown. 



The dependence of the transmission gap, A g , on the 
width of the ribbon with zigzag edges can be determined. 
For ribbons with armchair edges, in contrast, the gap is 
not well-defined because of the appearance of resonances 
and thus there is no clear indication of the gap's width. 
Shown in Fig. [5] is a plot of A g as a function of N for 
ribbons with zigzag edges. The length of the ribbons are 
m = 100. Also shown are the curves for linear, power-law 
and exponential fits. The values of the fitting parameters 
are shown in Table [TJ The correlation coefficients found 
for the three fits are close, with the exponential fit having 
the best value. 

Shown in Fig. [7] is the evolution of the wavefunction 
density \tp\ as the energy E is varied from within the 
transmission gap into the E>0 transmitting region of the 
N = 1, m = 8 ribbon with zigzag edges. In the figure, the 




Fig. 7. (color online) Wave function density in N = l, m = 8 
ribbons with zigzag edges. The energies are (a) E = 1.03, (b) 
£ = 1.05, (c) £ = 1.07, and (d) £ = 1.09. These energies are in 
the region around the right edge of the transmission gap. The 
colors aid the eye in determining the height of \tp\ 2 ■ 

ribbon is along the xy-plane. The input lead is coupled to 
the ribbon at x = and the output lead is coupled at the 
other side. For energies within the transmission gap the 
wave function is zero. For energies near the edge of the 
gap, for example at E — 1.03 as shown in Fig. [TJa), the 
wave function is almost zero everywhere except for a slight 
hump at the input side of the ribbon. This hump grows as 
the energy is increased until within the transmitting region 
a well-developed wave with a relatively large amplitude 
appears. Within the gap the wave is evanescent [35] and 
trying to penetrate the ribbon. As we increase the energy 
of the particle the wave eventually gains enough energy to 
proceed into the ribbon and transmit to the other side. 




Fig. 8. (color online) The |i/>| 2 in JV = 1, m = 8 ribbons with 
zigzag edges. These waves correspond to transmission peaks 
with energy values (a) £=1.12, (b) £ = 1.38, (c) £ = -1.52, 
and (d) £=-1.42. 

Shown in Fig. [8] are some of the wave function densi- 
ties that correspond to transmission peaks in iV = l, m = 8 
ribbons with zigzag edges. The wave in Fig. [5ta) is the 
highly transmitting wave that the waves in Fig. [7] develop 
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into as the particle's energy is increased. Compared to the 
waves in Fig. [7] the highly transmitting wave in Fig. Ufa) 
is sitting symmetrically at the center of the ribbon. This 
characteristic is a property satisfied by the wave functions 
that correspond to the transmission peaks in Fig. [31 in- 
cluding those for the armchair edges. 



4 Symmetry of ribbons with armchair edges 

In Fig. [3] the transmission coefficient for ribbons with arm- 
chair edges is shown to be symmetric with respect to a 
change in the sign of the energy. This symmetry arises 
from a property of the A matrix for the ribbon. 

Consider, for example, the linear problem Axp = b shown 
in Eq. (J4j) for the armchair hexagon. For a ribbon of N r 
sites the corresponding matrix A has size (N r + 2) x (N r + 2) 
and ip and b are column vectors of size (N r + 2). Suppose 
we have two linear problems, Atpi = b and Bif>2 = b, with 
the same right-hand side b. The and ip2 can be deter- 
mined using the inverses of A and B, i.e., = A _1 b and 
i/>2 = B _1 b. We can write the inverses of the matrices in 
terms of their corresponding matrix of cofactors, 



C T D T 
A andB=— , 



(5) 



where |A| and |B| are the determinants of the matrices. 
Since A and B are symmetric, C T = C and D T = D. In 
particular, can be calculated from 



V h J 
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Ja\ 



Cl,0 



CV r ,o 

\ CjVr+l, 



Co,N r +l \ 
Cl,W r +l 



CN r ,N r 



/2i/3\ 






V o / 



( 6 ) 

The transmission amplitude t\ can be determined from 
the multiplication of the last row of C with b. We get 



ti 



2i^Cjv r +i,o 



(7) 



A similar calculation can be done to determine t 2 . The 
transmission coefficients can then be calculated from T\ = 
|*i I 2 andT 2 = |i 2 | 2 . 

Suppose A=A(E) and B=A(-E). The condition that 
a change in the sign of E results in the same transmission 
coefficient, using Eq. ([7]), is 



|C(-E)jv r +i,o| — |C( _ E) Nr+ i >0 \ . 



(8) 



where we use the fact that |A(£:)| 2 = |A(-£0| 2 . For there 
to be symmetry between the positive and negative energy 
regions, Eq. (jSJ) must be satisfied. Note that this condition 
is independent of the size of the ribbon, as long as that 
ribbon is coupled to input and output leads. 

Calculating the cofactor C(E) 7 ^ of A in Eq. (|3|), noting 
that in our simulations m = v, results in 



This is an even function of E and thus satisfies the symme- 
try condition of Eq. ([5]). In contrast, for the zigzag hexagon 
its corresponding A matrix will have four more v entries 
than the one shown in Eq. ((4|). This reflects the second 
channel available in the coupling between the supersites 
and the hexagon. The cofactor of the resulting A for this 
case leads to 



|C(£) 7 ,o| 2 = Av 8 (E e 



4vE 5 + 2v 2 E 4 - 8v 3 E 3 
~7v 4 E 2 +4v 5 E + 4v 6 ). 



(10) 



\C{E) 



7,0| 



4v 10 (E 4 - 2v 2 E 2 



(9) 



Terms with odd powers of E appear and the cofactor 
therefore will be different when the sign of E is reversed. 
The symmetry condition in Eq. ||8]) is not satisfied and, 
as can be seen in Fig. [3fb) , the profile of the transmis- 
sion coefficient is different between the positive and nega- 
tive energy regions. This analysis can be extended to rib- 
bons with any N and m. As can be seen in Figs. 2] and [5] 
the symmetry is preserved only in ribbons with armchair 
edges. 



5 Summary 

We investigate the transport of a quantum particle travers- 
ing through honeycomb lattice ribbons with either arm- 
chair or zigzag edges with varying length and width. The 
ribbons are coupled to semi-infinite linear chains that serve 
as the input and output leads. We have seen the role of 
coupling the leads to the ribbons. In narrow, N — 1, rib- 
bons we find transmission gaps for both types of edges. 
In armchair edges the center of the gap coincides with 
the middle of the band. This symmetry is due to a prop- 
erty of the corresponding A matrices. In zigzag edges the 
gap center is displaced to the positive energy region. We 
also find transmission resonances for both types of edges. 
These transmission resonances occur whenever the wave 
fits symmetrically within the ribbon, reminiscent of reso- 
nant tunneling but without having finite potentials or tun- 
neling in this model. Increasing the length of the ribbons 
increases the number of transmission resonances without 
affecting the width of the transmission gap, as long as 
the ribbon is longer than a critical length when the gap 
can form. Increasing the width of the ribbons affects the 
width of the transmission gap. In ribbons with armchair 
edges transmission resonances appear within the gap as 
the ribbon width is extended. In ribbons with zigzag edges 
the gap systematically shrinks as the ribbon width is in- 
creased. 



We would like to thank Barbaros Ozyilmaz, Li- Fa Zhang and 
Jing-Hua Lan for insightful discussions. This work is supported 
in part by a Faculty Research Grant no. R-144-000-173-101/112 
from the National University of Singapore. 



References 

1. R. Saito, G. Dresselhaus, M.S. Dresselhaus, Physical Prop- 
erties of Carbon Nanotubes (Imperial College Press, Lon- 
don, 1998) 



E. Cuansing, J.-S. Wang: Quantum transport in honeycomb lattice ribbons with armchair and zigzag edges 7 



2. K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, 
Y. Zhang, S.V. Dubonos, I.V. Grigorieva, A. A. Firsov, Sci- 
ence 306, 666 (2004) 

3. K.S. Novoselov, D. Jiang, F. Schedin, T.J. Booth, V.V. 
Khotkevich, S.V. Morozov, A.K. Geim, P. Natl. Acad. Sci. 
USA 102, 10451 (2005) 

4. K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. 
Katsnelson, I.V. Grigorieva, S.V. Dubonos, A. A. Firsov, 
Nature (London) 438, 197 (2005) 

5. Y. Zhang, Y.W. Tan, H.L. Stormer, P. Kim, Nature (Lon- 
don) 438, 201 (2005) 

6. A.K. Geim, K.S. Novoselov, Nature Materials 6, 183 (2007) 

7. Y.W. Son, M.L. Cohen, S.G. Louie, Phys. Rev. Lett. 97, 
216803 (2006) 

8. M. Fujita, K. Wakabayashi, K. Nakada, K. Kusakabe, J. 
Phys. Soc. Jpn. 65, 1920 (1996) 

9. K. Nakada, M. Fujita, G. Dresselhaus, M.S. Dresselhaus, 
Phys. Rev. B 54, 17954 (1996) 

10. K. Wakabayashi, M. Fujita, H. Ajiki, M. Sigrist, Phys. Rev. 
B 59, 8271 (1999) 

11. M. Ezawa, Phys. Rev. B 73, 045432 (2006) 

12. L. Brey, H.A. Fertig, Phys. Rev. B 73, 235411 (2006) 

13. K.I. Sasaki, S. Murakami, R. Saito, J. Phys. Soc. Jpn. 75, 
074713 (2006) 

14. D.A. Abanin, P.A. Lee, L.S. Levitov, Phys. Rev. Lett. 96, 
176803 (2006) 

15. V. Barone, O. Hod, G.E. Scuseria, Nano Lett. 6, 2748 
(2006) 

16. A.H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. 
Novoselov, A.K. Geim, Rev. Mod. Phys. 81, 109 (2009) 

17. M.Y. Han, B. Ozyilmaz, Y. Zhang, P. Kim, Phys. Rev. 
Lett. 98, 206805 (2007) 

18. Z. Chen, Y.M. Lin, M.J. Rooks, P. Avouris, Physica E 40, 
228 (2007) 

19. C. Berger, Z. Song, X. Li, X. Wu, N. Brown, C. Naud, 
D. Mayou, T. Li, J. Hass, A. Marchenkov et al., Science 
312, 1191 (2006) 

20. X. Li, X. Wang, L. Zhang, S. Lee, H. Dai, Science 319, 
1229 (2008) 

21. B. Ozyilmaz, P. Jarillo-Herrero, D. Efetov, P. Kim, Appl. 
Phys. Lett 91, 192107 (2007) 

22. T. Otsubo, Y. Aso, K. Takimiya, Bull. Chem. Soc. Japan 
74, 1789 (2001) 

23. N. Robertson, C.A. McGowan, Chem. Soc. Rev. 32, 96 
(2003) 

24. R. Pearman, D. Bong, R. Breslow, GW. Flynn, Polym. 
Mater. Sci. Eng. 89, 213 (2003) 

25. D. Bong, I.W. Tarn, R. Breslow, J. Am. Chem. Soc. 126, 
11796 (2004) 

26. I.W. Tarn, J. Yan, R. Breslow, Org. Lett. 8, 183 (2006) 

27. D. Schultz, F. Biaso, A.R.M. Shahi, M. Geoffroy, K. Rissa- 
nen, L. Gagliardi, C.J. Cramer, J.R. Nitschke, Chem. Eur. 
J. 14, 7180 (2008) 

28. X. Guo, J.P Small, J.E. Klare, Y. Wang, M.S. Purewal, 
I.W. Tam, B.H. Hong, R. Caldwell, L. Huang, S. O'Brien 
et al., Science 311, 356 (2006) 

29. X. Guo, M. Myers, S. Xiao, M. Lefenfeld, R. Steiner, G.S. 
Tulevski, J. Tang, J. Baumert, F. Leibfarth, J.T. Yardley 
et al., P. Natl. Acad. Sci. USA 103, 11452 (2006) 

30. L. Venkataraman, J.E. Klare, I.W. Tam, C. Nuckolls, M.S. 
Hybertsen, M.L. Steigerwald, Nano Lett. 6, 458 (2006) 

31. Y.S. Park, A.C. Whalley, M. Kamenetska, M.L. Steiger- 
wald, M.S. Hybertsen, C. Nuckolls, L. Venkataraman, J. 
Am. Chem. Soc. 129, 15768 (2007) 



32. P.W. Anderson, Phys. Rev. 109, 1492 (1958) 

33. P.G. de Gennes, P. Lafore, J. Millot, J. Phys. Chem. Solids 
11, 105 (1959) 

34. S. Kirkpatrick, T.P Eggarter, Phys. Rev. B 6, 3598 (1972) 

35. P.F. Bagwell, Phys. Rev. B 41, 10354 (1990) 

36. H. Schomerus, Phys. Rev. B 76, 045433 (2007) 

37. Y.M. Blanter, I. Martin, Phys. Rev. B 76, 155433 (2007) 

38. J.P. Robinson, H. Schomerus, Phys. Rev. B 76, 115430 
(2007) 

39. D. Daboul, I. Chang, A. Aharony, Eur. Phys. J. B 16, 303 
(2000) 

40. J.S. Wang, J. Wang, J.T. Lii, Eur. Phys. J. B 62, 381 
(2008) 

41. G.H. Golub, C.F. van Loan, Matrix Computations, 
3rd edn. (Johns Hopkins Univ. Press, Baltimore, 1996) 

42. E. Cuansing, H. Nakanishi, Physica A 387, 806 (2008) 

43. M.F. Islam, H. Nakanishi, Phys. Rev. E 77, 061109 (2008) 

44. E. Cuansing, H. Nakanishi, Phys. Rev. E 70, 066142 (2004) 



